Signature of antiferromagnetic long-range order in the optical spectrum 
of strongly correlated electron systems 
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We show how the onset of a non-Slater antiferromagnetic ordering in a correlated material can 
be detected by optical spectroscopy. Using dynamical mean-field theory we identify two distinctive 
features: The antiferromagnetic ordering is associated with an enhanced spectral weight above the 
optical gap, and well separated spin-polaron peaks emerge in the optical spectrum. Both features 
are indeed observed in LaSrMn0 4 [Gossling et al, Phys. Rev. B 77, 035109 (2008)]. 

PACS numbers: 71.27. +a, 71.10.Fd 
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I. INTRODUCTION 

This article deals with the following issue: Is there 
a distinctive signature of the onset of antiferromag- 
netic (AF) long-range order in the optical spectrum of 
a strongly correlated insulating antiferromagnet, for ex- 
ample an antiferromagnetic transition-metal oxide? 

AF ordering is a characteristic feature of strongly cor- 
related Mott insulators where superexchange drives the 
ordering of localized magnetic moments. On the other 
hand, AF can also arise from Fermi-surface nesting in 
a weakly correlated material. From a theoretical point 
of view, a continuous crossover driven by the interaction 
strength connects weak- and strong-coupling antifcrro- 
magnets, and an unambiguous distinction between local- 
moments and Fermi-surface AF is lacking. 

In this work we propose that, remarkably, optical spec- 
troscopy can be used to infer the correlated nature of 
an AF state. Optical spectroscopy is an invaluable ex- 
perimental probe of correlated materials, providing key 
physical information such as the optical gap, the rela- 
tive weight of low-energy Drude excitations in metallic 
systems, and quite importantly the transfers of spectral 
weights often observed in correlated materials when tem- 
perature or composition are varied (for a recent review 
see Ref. 0). Spin de grecs of freedom are however much 
less coupled to light than the charge so that the signa- 
tures of magnetism are expected to be comparatively 
weak. Zone-center magnons, for example, do not show 
up in the optical conductivity unless inversion symme- 
try is broken 1 -. Furthermore, the key energy scale as- 
sociated with magnetic ordering, the superexchange J 
(~ D 2 /U with D being half the bandwidth), is much 
smaller than the scale U corresponding to the local ma- 
trix element of the screened Coulomb interaction, and 
hence generally significantly smaller than the optical gap 
itself (~ U - 2D). 



Despite these shortcomings, we show that clear signa- 
tures associated with antifcrromagnetism are expected 
in the optical conductivity below the Neel temperature 
(Tjv) when the system is in the strong-coupling (Mott 
insulating) regime. These signatures are twofold. First, 
there is a spectral weight transfer as one cools the system 
from above the Neel temperature down to low tempera- 
ture. This transfer actually provides a diagnostics of the 
strong-coupling (superexchange) or weak-coupling (spin- 
density wave) nature of the antiferromagnetism. In the 
former case, spectral weight increase is observed, corre- 
sponding to a kinetic energy gain associated with the 
ordering. In the latter case, a spectral weight decrease 
takes place, corresponding to a kinetic energy loss (poten- 
tial energy gain). For superconducting long-range order, 
similar diagnostics of weak vs. strong-coupling pairing 
have been discussed and probed experimentally in the 
context of cuprate superconductors 2 and for models with 
attractive interaction 3 -^. Second, a multi-peak structure 
develops above the gap in the ordered phase at strong 
coupling. These peaks correspond to the spin-polaron 
excitations associated with the motion of a hole in the an- 
tiferromagnetic background. They have been previously 
discussed in the context of on&jjarticle spectroscopy such 
as photoemission in Refs. (aQ and this work identifies 
their signature in the optical conductivity. 

Finally, we consider the possibility of detecting these 
signatures experimentally, and discuss in detail the case 
of LaSrMn04 where we argue that some of these effects 
may already have been observed^. 

The outline of the paper is as follows: In Section |TT] 
we summarize the main aspects of our calculation of the 
optical conductivity for the Hubbard and t-J model. De- 
tails on the DMFT calculations for the two models are 
provided in Appendix I VI Al and IVIB1 respectively. In 
Section HTT1 we discuss the main theoretical results before 
we compare them to experiment in Section IIVI A sum- 
mary and conclusion is given in Section El 
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II. MODELS AND METHODS 

In order to study the optical conductivity of strongly 
correlated antiferromagnets, we perform dynamical 
mean-field theory (DMFT) calculations for a half-filled 
Hubbard model, described by the following Hamiltonian: 

H = -t J2 4 c i- + U Y1 "W- C 1 ) 

(ij)a i 

Here, t denotes the hopping amplitude between nearest- 
neighbors, U the Coulomb interaction, and c\ a {ci a ) cre- 
ates (annihilates) an electron with spin a on site i; 
v,i a = c| cr Cj CT - For simplicity, we use a semi-elliptical den- 
sity of states of half-width D, corresponding to a Bethe 
lattice with infinite coordination. Our results however 
hardly depend on this specific choice. We vary the tem- 
perature T from above to below the Neel temperature. 
We therefore need the paramagnetic and the antiferro- 



magnetically ordered solutions, which can both be calcu- 
lated within DMFTM. 

It is known that the Hubbard model for large values 
of the local repulsion U maps onto the t — J model. For 
the specific model under consideration, the relation be- 
tween J and U is the following: for an individual pair of 
sites we have J- ln d. = 4t 2 /U and summed over all nearest 
neighbors z we have J = zJ- m< \. = D 2 /U (coined J* in 
Ref. d) . The quantities of interest in the t—J model in in- 
finite dimensions can be expressed in the form of a simple 
continued-fraction^^ while the Hubbard model requires 
a more computationally expensive solver. Therefore we 
will often use the t—J solutions in the case of the antifcr- 
romagnetic phase. Indeed, the electron-removal part of 
the fc-integrated spectral function of an half-filled Hub- 
bard model below Tjy coincides up to higher order terms 
in J with the spectral function of a single hole in the t—J 
model. The way we compute the optical conductivity for 
the t — J model, and in particular how we get across the 
absence of an electron addition part of the spectrum of 
the t — J model, is explained in detail in the Appendix. 



Within DMFT the optical conductivity can be expressed exactly in terms of the convolution of two single-particle 
Green functions without vertex correction o 9 ^ 1 1 1 12 . In the antiferromagnetic phase this reads^: 
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[lmGf(e,w)ImGf(e,w + SI) + ImGf (e,w)ImGf (e,u + SI)] , 



(2) 



where / denotes the Fermi distribution function, G AA and G AB are respectively the Green function of an electron 
with spin a hopping between sites of the same sublattice and between sites of the two different sublattices A and B, 
in which the Bethe lattice is split in the AF phase. 

I 



In this expression, the quantity -^^(D 2 — e 2 ) 3 / 2 repre- 
sents the product of the density of states of the Bethe 
lattice times the current vertex. This assumption insures 
that the f-sum rule associated with the optical conduc- 
tivity holds, namely that its integral is proportional to 
the kinetic energ y 11 ! 14 " — . 



III. MAIN RESULTS 

In Fig. [T] we show our main result, namely the 
temperature-dependence of the optical conductivity as 
one cools the system from close to the Neel tempera- 
ture down to low temperature, in a regime of strong 
correlations, i.e., U/D = 5. For temperatures slightly 
lower than Tjy, the optical conductivity resembles that 
of the paramagnetic insulator: it has a gap of order 
U — 2D, and the absorption peak above the gap is essen- 
tially featureless. Upon cooling, two remarkable features 
characterize the optical signal: (i) the optical absorp- 
tion strongly increases and (ii) a multi-peak structure 
emerges. Such a structure becomes more and more vis- 
ible as one goes down in temperature. Note that these 



changes occur without a visible variation in the position 
of the absorption edge (size of the gap). As discussed be- 
low, they should nevertheless be detectable by accurate 
experiments as a function of temperature. 

Later in this section we will explain the physical origin 
of the peaks in the optical conductivity by investigating 
the single particle spectral function and its evolution with 
U. At this stage, let us only stress that these peaks are 
by no means an artifact of the discretization required by 
our exact-diagonalization DMFT solver—. To better com- 
pare with experiments, we used a Lorenztian broadening 
of 0.05-D to plot the optical conductivity. This value is 
small enough to distinguish the first peaks but also large 
enough to mimic the smearing of the multi-peak struc- 
ture due to experimental resolution. 

We now discuss in more details the two key effects we 
mentioned before. 

i) Change of spectral weight through Tjy. The increase 
of spectral weight upon cooling is characteristic of a Mott 
antiferromagnet at strong coupling. Therefore it is the 
first gross feature one may look for in an optical exper- 
iment in order to establish whether the antiferromag- 
netism in the material is of strong-coupling nature, with 
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FIG. 1. (Color online) Temperature evolution of the optical 
conductivity of the t — J model at strong coupling (U = 5D). 
Inset: temperature dependence of the absolute value of the 
kinetic energy (— K in units of D) and of the staggered AF 
magnetization m. 



preformed magnetic moments, or it is of weak-coupling 
nature, with a collective Fermi-surface instability. To un- 
derstand this, we can use the fact that the integral of the 
spectral weight is directly proportional to the kinetic en- 
ergy. At strong coupling U 3> D, the ordered phase is 
stabilized by a gain in kinetic energy, associated to the 
establishment of long-range coherence of the local mag- 
netic moments. In simple words the system gains kinetic 
energy by increasing the staggered magnetization: Co- 
herent hopping processes take place in a staggered spin 
background are favored if the spin pattern is as close as 
possible to the Neel one 3 -^. Thus the lower the tempera- 
ture, the larger the staggered magnetization and in turn 
the larger the spectral weight. 

If we consider instead the weak-coupling regime U <C 
D, a totally different mechanism for the stabilization of 
the ordered phase takes place. The AF phase is stabi- 
lized by a (small) potential energy gain, as the onset of 
a non-zero order parameter is correctly described by the 
static mean field theory in this regime. These theoretical 
considerations are well confirmed by our numerical data 
shown in the bottom panel of Fig. [5] where we compare 
K AF and K PM , the kinetic energy of the ordered and 
disordered phase respectively. The complete evolution 
of the kinetic and potential energy as a function of U is 
shown in Fig. [3J At strong coupling the antiferromag- 
netic order is clearly stabilized by the kinetic energy while 
at weak coupling (U < 2D) the ordered phase is driven 
by the gain in potential energy. At intermediate values 
of the interaction strength, corresponding to the cross- 
over between weak and strong coupling, we have a mixed 
regime in which both energy differences are negative 3 -^, 
i.e. where the onset of the AF is stabilized by both ki- 
netic and potential energy. 

ii) Spin-polaron peaks. To clarify the origin of the 
multi-peak structure, we compare in Fig. [5] the evolu- 
tion from weak to strong coupling of the optical con- 



ductivity to that of the one-particle spectral function. 
It is evident that the multi-peak structure of the spec- 
tral function^ survives in the optical conductivity. The 
peculiar behavior of the spectral functions was first un- 
derstood in a pioneering work by Strack and Vollhardt 5 
for the t—J model in infinite dimensions; and by some of 
us more recently within the Hubbard model in the anti- 
ferromagnetic phase by means of DMFT 6 . The multiple 
peaks arise from string-like excitations associated with 
a hole moving in a Neel background. In infinite dimen- 
sions the lifetime of such string excitations is infinite due 
to the absence of quantum spin-flip processe a 5 ' 6 ' 19 which 
would act as 'damage-repairing' processes, which means 
that the Heisenberg interaction reduces to the Ising one 
(see Appendix IVI Al for further details). This results 
in the formation of coherent but heavy quasiparticles, 
"dressed" by a spinon cloud which is commonly called a 
"spin-polaron" . 

The formation of the spin-polaron can be described 
only using methods (such as DMFT) that go beyond the 
static mean field of a spin-density wave in the mean- 
field Hartree-Fock approximation. This is illustrated in 
the upper panels of Fig. [5J where we compare the evolu- 
tion with U of the dynamical (red line) and static mean- 
field (green line) results for the electron-removal spec- 
trum. One clearly sees that the only feature of the static 
mean-field (Hartree-Fock) calculation is a Slater-like sin- 
gle peak centered around —{7/2. The peak has weight 
Z = 1 and its width is given by </, i.e. the width shrinks 
to as 1/ U upon increasing U. In the dynamical mean- 
field calculation, instead, the weight of the lowest lying 
excitation is much smaller (Z is proportional to J/D) 
and the remaining weight is transferred to the higher- 
energy peaks. This makes a total bandwidth of order 2D 
rather then J, resulting in a smaller gap compared to 
the Hartree-Fock calculations. 6 . This is why the gap in 
Fig. [1] is of order U — 2D, while it would be exactly U 
within static mean-field theory. 

In the t — J model in infinite dimensions (as well as 
in the Hubbard model within AF-DMFT, up to higher- 
order corrections in t 2 /U) the position of the peaks can 
be analytically expressed through the zeros of the Airy 
functions^ and the spacing between them goes as J 2 / 3 . 
As we will also discuss in Section [V] these analytic ex- 
pressions hold strictly only in infinite dimensions. Hence 
the applicability of our results to real materials, and in 
particular the emergence of spin-polaron peaks, strongly 
relies on the assumption that it is possible to neglect the 
quantum spin-flip processes. In principle this is not guar- 
anteed for finite dimensional systems, therefore one may 
expect some of the peaks to broaden or even to be washed 
out, due to the coupling to dispersive spin wave s 6 ' 22 . 
On the other hand, recent sophisticated calculations for 
the two-dimensional t — J model give clear evidence for 
the survival of the first one or two spin-polaron peaks 
with a similar separation as in the infinite-dimensional 
calculation 20 . Furthermore the case of a double layer an- 
tiferromagnet was recently considered in Ref. [2ll There, 
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FIG. 2. (Color online) Lower panel: Evolution of the optical conductivity of Hubbard model from weak (left) to strong 
coupling (right). Shown are both, AF phase and PM phase. In the figure the corresponding values of the optical integrals 
(=kinetic energy) are also reported, showing a change of the "hierarchy" between the kinetic energy scales of the two phases 
at intermediate coupling. 

Upper panel: Corresponding spectral function with a multi-peak (spin-polaron) structure.— 
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FIG. 3. (Color online) Difference of kinetic, potential and 
total energy between AF and PM phase from weak to strong 
coupling. At strong coupling AF is stabilized by the kinetic 
energy and at weak coupling by the potential energy. 



it is shown that the inter-layer exciton made of one hole 
and one doublon in each layer displays confinement ef- 
fects when the inter-layer coupling is smaller than the 
intra-layer one. In such a situation, the spin-polaron exci- 
tations discussed here may become particularly relevant. 

When considering three dimensional cases, where 
DMFT is typically more accurate and finite- 
dimensionality effects are definitively weaker than 
in two dimensions 23 , spin-polaron features are expected 



to be even more visible. In the following section, we 
discuss indeed the relevance of our results for the three 
dimensional manganite LaSrMnC>4. 



IV. COMPARISON WITH EXPERIMENTS 

On the basis of the previous discussion, the ideal ma- 
terial for the observation of the temperature evolution 
predicted above would be a completely isotropic three- 
dimensional Mott antiferromagnet. To our knowledge 
the most relevant set of measurements in this respect 
is that by Gossling, et al.J on LaSrMn04, a three- 
dimensional C-type layered antiferromagnetic material 
with Tjv = 133K, while we are not aware of similar stud- 
ies as a function of temperature on G-type antiferromag- 
nets, which would be the closest realization of our large- 
coordination results. In Fig. |4]we show the data of Rcf. 
7 for the optical conductivity of LaSrMnCU along the a 
axis at different temperatures. The two main structures 
around 3.5 eV and 4.5 eV have been explained in terms 
of multiplet split transitions from a d A to a d 5 configura- 
tion of the Mn. Our main interest instead is focused on 
the fine structure arising upon cooling on top of both of 
these main structures. 

In the light of our theoretical analysis, we can di- 
rectly relate this observation with the physics of the spin- 
polarons described in the previous section. Since the two 
structures are well separated, we shall make a simplify- 
ing assumption (not meant to be quantitative, but suf- 
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FIG. 4. (Color online) Optical conductivity of LaSrMnCU 
between 0.75 and 5.80 eV for different temperatures. The 
raw data have been kindly provided by the authors of Ref. 7. 

ficient to reveal the main qualitative aspect of the phe- 
nomenon). Namely, we shall assume that we can adopt 
for each d A —> d 5 transition a simple single-band Hubbard 
model description. 

Since the separation between the multiplet peaks is 
larger than the energy scale associated with the spin- 
polaron excitations, we expect to be able to detect the 
above-mentioned increase of spectral weight around each 
of the main absorption peaks and the consequent ap- 
pearance of spin-polaron structures when going below the 
Neel temperature. 

All curves below TV display, in the same region, a very 
similar low-frequency tail (below ~ 3 eV) while in that re- 
gion the higher-temperature measurements show a grad- 
ual increase of spectral weight that we attribute to a more 
standard thermal broadening effect. It is likely that such 
thermal effects are visible also in the low-energy multi- 
plet peak at 3.5 eV, which lies closer to this tail. This 
makes this spectral feature less visible with a first dis- 
tinguishable excitation shifting with temperature unlike 
what happens for the second multiplet structure at 4.5 
eV. 

We focus therefore on the high-energy spectral struc- 
ture around 4.5 eV, which seems to be unaffected by the 
temperature broadening of the low-frequency data, and 
we zoom the plot to focus on this feature. In order to 
highlight the temperature effect, we display in Fig. [5]the 
difference between the optical conductivity at a given 
temperature and the conductivity at the highest mea- 
sured temperature (300K). 

Fig. [5] and its theoretical counterpart (Fig. [6]) en- 
dorse our physical interpretation of the peaks in the op- 
tical signal. Although the frequency grid in the experi- 
mental data is 0.05 eV, at least two well-spaced peak- 
structures emerge below TV and, as expected within 
our spin-polaron picture, the lower the temperature the 
weaker the temperature dependence becomes. Indeed be- 
tween 52 and 15K there is hardly any difference in the 
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FIG. 5. (Color online) Difference between the optical con- 
ductivity of LaSrMn04 between 4.0 and 5.5 eV at differ- 
ent temperatures with respect to the optical conductivity at 
T = 300K. Data by Gossling, et al.Z 
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FIG. 6. (Color online) Difference between the computed op- 
tical conductivity at different temperatures and the one at 
TV,— showing a very similar structure as the experiment in 
Fig. El 



optical conductivity in this frequency region, as one can 
expect from the behavior of the staggered magnetization 
(see inset to Fig. [I}, which is almost saturated far below 
T/v- The appearance of the spin-polaron fine structures 
become clear below Tn = 133K. Yet, the polaronic fea- 
tures start to be visible already in the 150K curve. We at- 
tribute this to fluctuation effects above TV: even if long- 
range order is not yet fully developed, a sufficiently long 
coherence length of the spin-spin correlations is enough 
to sustain a non-fully developed spin-polaron mechanism. 

Similar features as the ones we have described are 
present in the optical absorption of another material: 
KCUF3. For temperatures lower than the magnetic or- 
dering temperature of this compound a multi-peak struc- 
ture is indeed observed (see Fig. 2b of Ref. 24), for which 
an interpretation in terms of spin-waves has been given. 



6 



We have therefore found in nature what we expect from 
our model calculation: A clear signature of antifcrromag- 
netism arising from strong correlations in the optical con- 
ductivity of a strongly correlated material which emerges 
as the system is cooled below Tn. If the material is suffi- 
ciently three-dimensional and if the resolution of the ex- 
periment is good enough, a clear signal with a multi-peak 
structure can be observed. Another way of revealing the 
same effect would be to compare the increase of spec- 
tral weight to the measured staggered magnetization as 
a function of T. 



V. CONCLUSIONS 

We have calculated the optical conductivity of a 
strongly correlated antiferromagnet described by the 
Hubbard or the t — J model. Within dynamical mean- 
field theory, we have shown that the optical conductivity 
is characterized by the presence of equally spaced spin- 
polaron peaks which appear right below the Neel tem- 
perature and become more and more pronounced as the 
temperature is reduced. The onset of antiferromagnetic 
ordering is stabilized by a kinetic energy gain in strong 
coupling superexchange magnets, while at weak-coupling 
antiferromagnetism is due to Fermi-surface nesting and 
it is stabilized by a static mean-field like potential energy 
gain. 

In the strong coupling regime, the sharp peaks associ- 
ated to spin-polarons originate from string-like excitation 
created by the motion of a hole in an antiferromagnetic 
background, and we can safely distinguish them from ar- 
tificially sharp peaks characteristic of exact diagonaliza- 
tion solutions of DMFT. 

Even if our formal derivation of the spin polaron peak 
relies on the infinite-coordination limit in which DMFT 
becomes exact, we expect at least part of these features 
to survive in finite dimensions, and in particular in the 
three dimensional case, in which DMFT typically proves 
accurate^. We notice that at least the lowest-energy ex- 
citation peak is found also in two dimensions, as shown by 
different approaches--, testifying that the physics behind 
our DMFT results survives when the large dimensionality 
limit is far. 

If we look for experimental realizations of this scenario, 
we should take into account temperature effects and the 
finite experimental resolution, which can further smear 
out some of the fine structures. In this light, it is remark- 
able that an experimental study on LaMnSrO^ reveals 
a multi-peak structure whose form and temperature be- 
havior are in very good agreement with our calculations 
on top of the multiplet structure. 

Moreover, considering the recent work by Rademaker 
et al£^, optical spectroscopy performed along the c-axis 
in layered systems may be a promising tool to detect and 
investigate the occurrence of the spin-polaron physics. 

An even stronger signatures of the spin-polaron effects 
are to be expected in three dimensional correlated ma- 



terials which present a three-dimensional G-type antifer- 
romagnetism, like SrMnC>3, CaMnC>3, BaMn03, calling 
for future IR-spectroscopy investigations of the PM-AF 
transitions in these compounds. 
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VI. APPENDIX 

In this Appendix we give all details about how to con- 
struct an approximated expression for the DMFT self- 
energy of the half-filled Hubbard model on the basis of 
the DMFT solution for the corresponding t— J model both 
in the antiferromagnetic and the paramagnetic phase. 
We will discuss also the validity of this procedure by ex- 
ploiting the mapping of the (half-filled) Hubbard Hamil- 
tonian onto the t — J one2£ for U ~S> D and by comparing 
our results for the self-energy and the optical spectra ex- 
tracted from the DMFT solution of the t — J model with 
the correspondent ones for the Hubbard model. The lat- 
ter are available with high enough numerical precision. 



A. DMFT for the t — J model 

The main motivation to use the t — J results for ap- 
proximating the DMFT self-energy stems, as mentioned 
in Sec. II, from the simplicity of the algorithm for solving 
the DMFT problem of one hole in the t— J model with 
(complete or partial) AF-Neel order— as opposed to the 
bigger numerical effort needed for the Hubbard model. 
Specifically for one-hole at T = one finds (for more 
details see Ref. [5|) that for d — > oo 

1. the quantum spin-fluctuations associated to the 
term proportional to SfSf + SfSlj of the SU(2)- 
symmctric Heisenberg interaction are completely 
frozen. In other words, the Heisenberg coupling 
of the t—J Hamiltonian is reduced to an Ising cou- 
pling, whose ground state is the AF-Neel state with 
a perfectly staggered spin configuration, i.e, with 
full staggered magnetization (m = J2 r . ( n it — 
n i,i){~l) Ti — lj where N represents the total num- 
ber of sites); 
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2. the "retraceable path" approximation becomes rig- 
orously exact: In a perfect AF-Neel state, the local 
physics (e.g., the local Green function G(ui)) is built 
only by hopping processes in which the one hole 
comes back to the starting site following precisely 
the same path. In fact, any non-retraceable path 
is associated with an higher-energy cost, as the 
hole movement would result in a path of misaligned 
spins in the AF-Neel background. This implies that 
the geometry of the underlying lattice becomes ir- 
relevant in d — > oo: If only retraceable paths mat- 
ter, every geometry becomes topologically equiva- 
lent to the Bethe lattice, a tree-structure which by 
construction only allows for retraceable paths. 

As a consequence at T — the discrete energy lev- 
els characterizing the analysis of one-hole dynamics in 
the AF-Neel background of the t — J model can be re- 
lated to the simple quantum mechanical problem of one 
particle in a step potential V(n) = nJ, being n the 
number of misaligned spin created by the hopping pro- 
cesses of the hole. This allows in turn to obtain a com- 
pact continued-fraction expression for the retarded Green 
function G t -j(uj) = G^j(oj + i0+) (see Ref. H) 



of P, when smaller than 1, will affect the self-consistent 
expression for the retarded Green function: 



Gt-J(LO) 



D 2 1 ' 

- — G w (w - - J) 

1 

"752 i 



(3) 
(4) 



which can be solved with high numerical accuracy by 
simply truncating the continued fraction to the desired 
level of precision, and then building up the full expres- 
sion iteratively. Complications arise at finite T: In par- 
ticular point 2 is in principle no longer valid, as thermal 
fluctuations will reduce the value of the staggered mag- 
netization, allowing for a contribution of non-retraceable 
paths to the local physics. However, we can easily cir- 
cumvent this problem, if we consider explicitly the Bethe 
lattice case for our DMFT calculations. This is, on the 
one hand, a reasonable choice, as the finite bandwidth 
of the Bethe lattice reflects more directly the realistic 
situation of finite bandwidths in real materials, and, on 
the other hand, it preserves by construction the valid- 
ity of the "retraceable path" condition at all tempera- 
tures. Working with the Bethe lattice allows, in turn, for 
a straightforward generalization of the continued-fraction 
expression for G(ui) to the finite temperature case, along 
the lines of the exact derivation of ReffTol. Specifically, 
one can easily generalize Eq. ([3]) to finite temperatures 
by taking into account the probability P (or 1 — P) of 
finding correctly aligned (or misaligned) spins along the 
retraceable path of the hole. The probability P as a 
function of T is easily expressed as P = |(1 + m(T)) in 
terms of the staggered magnetization m(T). The latter, 
because of point 1 above, can be directly obtained by 
the simple mean field Curie- Weiss self-consistent equa- 
tion m(T) = tanh(^m), with T N = \j . The value 



D 2 



uj +— (PG f _,/(w-wp) 



+ (l-P)G t ^j(u + io P ))}- 1 



(5) 



being up = up(T) = ^m(T)J, and ui = uj + i0 + . Note 
that this expression is exactly reduced to Eq. (|3j) in the 
limit T ->■ 0, as ni(T = 0) = P(T = 0) = I. The 
additional terms in Eq. ([3]) at finite T result in a in- 
creasing proliferation of new (smaller) peaks no longer 
exactly matching with the original multi-peak structure 
of the continued- fraction expression at T = 0. As a con- 
sequence, by increasing T, one observes a gradual smear- 
ing out of the S-\ike peak structure of Eq. ([3]), which 
is completely washed out at T > T/v, when one has 
m = 0, P = 5 and the spectral representation G(oj) is re- 
duced to the characteristic semicircular shape of the non- 
interacting Bethe-lattice case. This result can be eas- 
ily understood, because the Heisenberg interaction term, 
which is non-local in the site indices becomes irrele- 
vant due to the d — > oo scaling in the paramagnetic phase 
(i.e., for m = 0) and, therefore, for T > TV one gets the 
same spectral function of the case J — 0, independently 
from the actual value of J. 



B. Approximation for the DMFT self-energy of the 
Hubbard model 



After having recalled how to solve exactly the DMFT 
problem of one hole in the t—J model, we will discuss here 
how to approximate the DMFT self-energy for the Hub- 
bard model in the limit of intermediate-to-strong cou- 
pling by exploiting this solution. In fact, it is known^ 
that in the limit U D the Hubbard model can be 
mapped onto a corresponding t — J model with J = — . 
When using the t — J results to approximate the Hub- 
bard DMFT self-energy or the k— integrated spectral and 
Green functions in the limit of U 3> D, one should not 
forget that the exact "mapping" between the two models 
is obtained by projecting out the doubly occupied sites. 
This evidently means that the t — J model cannot de- 
scribe the whole frequency range of the Hubbard model, 
even in the limit of U 3> D, but, at most, only "half" 
of it. This becomes quite obvious, when explicitly con- 
sidering the k— integrated spectral functions: The gap 
between the Hubbard bands is originated by the energy 
difference (~ U) between single occupied and double oc- 
cupied/empty sites, but this energy difference is exactly 
what is projected out in the mapping onto the t—J. 
Therefore, when computing the DMFT Green function 
for one hole in the t — J model, the associated spectral 
function will be related to electron-removal (uj < 0, i.e. 
photoemission) part of the k— integrated spectral func- 
tion of the Hubbard model: More precisely, taking into 
account the above mentioned energy shift of order U (not 



8 



0.6 



0.2 
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FIG. 7. (Color online) Comparison between the staggered 
magnetizations m as a function of temperature obtained 
by means DMFT calculations for the Hubbard model (with 
U — 5D) and for the corresponding t — J model: Due to the 
residual (but still finite) presence of double occupied/empty 
sites in the AF-ground state of the Hubbard model, the mag- 
netization (as well as the value of Tjv) is slightly lower in the 
former case, and -at difference with the t—J case- it never 
reaches 1 (complete magnetization), even for T — > 0. 



explicitly included in the t — J Hamiltonian), it will cor- 
respond to the lower Hubbard band. Similarly, if one 
considers the problem of one electron added to the t — J 
model, the associated spectral function would correspond 
to the upper Hubbard band of the Hubbard model spec- 
trum. 

These considerations are easily translated in a simple 
practical implementation: In fact, considering the sym- 
metry with respect to the change u> <H> — uj for both the 
paramagnetic and the antiferromagnetic case, one can 
obtain the t— J- like approximation of the k— integrated 
retarded Hubbard Green function G a (io + i0 + ) on a given 
sublattice (for the sake of defmiteness we identify with 
a =t the majority spin component on the AF order on 
the chosen sublattice) directly from the retarded Green 
function (or the spectrum) Gt-j(ui) (A t -j(u)) of the one- 
hole t — J problem as 

G t (oj)= PGW(w+|) + (l-P)Gw(^+j) (6) 

G^(«) = (l-P)G! w (a;+|)+PGfw(-w+|) (7) 

The essential ingredients of Eqs. ([BJ and (0 are, hence, 
two: (i) the spectral (Green) function of the one-hole 
problem of the t — J model, and (ii) the probability P to 
find (on the chosen site) the spin correctly aligned with 
the AF-Neel background. A word of caution is due at 
this point: As the value of P is defined in terms of the 
staggered magnetization m(T), the latter naturally rep- 
resents the most appropriate "link" between the DMFT 
t — J and Hubbard Green functions. In fact, in using the 
t — J approximation for the DMFT spectra of Hubbard 
model, one should keep in mind that the mapping is rig- 
orously exact only in the limit of U — > oo, otherwise, cor- 
rections ~ jj? have to be expected: From a more physical 
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FIG. 8. (Color online) Temperature evolution of the DMFT 
self-energy of the Hubbard model, calculated from the DMFT 
Green function of the t — J model, according to Eqs. (JTJl and 
|9]for the specific case U — 5D(J — 0.2D). Shown is the spin- 
up self energy on the sites with spin-up majority; the spin- 
down self energy is the particle-hole counterpart (obtained 
by mirroring uj — s> -co). As for the requirement of an equal 
staggered magnetization between the Hubbard and the t—J 
model, the corresponding temperature of the t — J case are 
reported in the legend as Tt-j . The dotted vertical lines are 
located at uj = -|m and uj — +^rn i.e., where in a static 
mean-field context one would find the two peaks of the self 
energy. 



perspective, the presence of doubly occupied/empty sites 
in the ground state of the Hubbard model stays indeed 
finite (though becoming smaller and smaller for increas- 
ing U) at any value of the Hubbard interaction. As a 
consequence, for any given set of parameters U and T, 
the staggered magnetization of the Hubbard model will 
be slightly lower than that for the t — J model, as it is 
also found in our DMFT data, plotted in Fig. [7J 

It should be clear, hence, that the best agreement be- 
tween the Hubbard model (for a given U) and the corre- 
sponding t—J (with J = -jj-) approximation is found for 
sets of Hubbard and t — J Green/spectral function with 
the same magnetization values (which are achieved, ac- 
cording to Fig. [71 for slightly different temperatures^). 

Having defined the equations (and the conditions) for 
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best approximating the Hubbard DMFT Green functions 
with the t — J model ones, it is straightforward to ex- 
tract from these the correspondent DMFT expressions 
for the (spin-dependent) Hubbard self-energies £ CT , via 
the DMFT self-consistency equations (generally written 
to include the AF order 9 ): 



D 2 

S t (w) = oj + h- — [G±(u>)] 1 
D 2 

E 4 (w) = uj + ii- — [G t {uj)} \ 



(8) 
(9) 



where the chemical potential /i is set fixed to ^ (particle- 
hole symmetric case). Representative results of our 
DMFT t — J approximation for the imaginary part of 
our Hubbard self-energy ImTi^(uj), computed according 
to Eqs. ((6|) - (j9|) are reported in Fig. [8] The main fea- 
ture of the T— dependence of Itl\Yi^{uj) well agree with 
the corresponding DMFT data for the Hubbard model 
shown in Ref . 0, supporting the validity of the t — J ap- 
proximation. In particular, if U 3> D, both in Hubbard, 
as well as in its t — J approximation, one finds that for 
T <C T/v (upper panel of Fig. |5} ImX^(oj) is character- 
ized by a multi-peak structure for uj < which reflects 
that of the spin-up spectral function, and, hence, the 
above-mentioned spin-polaron physics. Remarkably the 
"magnitude" of ImK(w), e.g. if quantified by the value 
of its frequency integral, is not particularly large for low- 
T, implying a certain degree of "coherence" of the spin 
polaron excitations in this regime^. Note that, except 
for the case of full polarization a much weaker multi-peak 
structure is also visible for u) > 0, roughly speaking where 
the upper Hubbard would be located in the PM phase. 

By increasing T, one observes a rapid strengthening of 
such "secondary" multi-peak structure at u > 0, which 
becomes quickly predominantly (second and third panel 
of Fig. |5]). More specifically, one notes that the lower 
energy peak of the "secondary" (uj > 0) structure of 
ImE(tj), is strongly enhanced with increasing T, and at 
the same time undergoes a constant "softening", as its 
position (inside the spectral gap) tends to go to zero fre- 
quency for T — > Tn (lower panel of Fig. |5J). In the limit 
of T = T/v, when Y>^(lu) = £|(w), this (now extremely 
strong peak) would be located exactly in the center (that 
is at uj = 0), ensuring the existence of the Mott- Hubbard 
gap ~ U — 2D also in the absence of any long range 
magnetic order. This can be compared with the exact 
solution of the atomic limit (D = 0), for which one can 
easily compute the Green function^, e.g. for the spin t 
sector 
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and, from this (and the corresponding Gi(lu) — G^(— uj)), 
calculate the self-energy 
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FIG. 9. (Color online) Comparison between the DMFT opti- 
cal conductivity of Hubbard model (impurity solver: Lanczos 
with N s — 25) and corresponding t— J model. 



which evidently displays a peak in its imaginary part lo- 



cated at 



UJpeak = m- 



(i.e, a "softening" for increasing 



T/decreasing m) of weight w pea k — 7r^4-[l — m 2 ] (i.e., of 
increasing strength for increasing T/decreasing to). At 

T = T/v, one gets eventually Ef (a;) = Ej.(w) — — — 



4(w+i0+) ' 

whose divergence in uj = directly corresponds to the 
spectral gap of (exact) size U. In comparison with the 
atomic limit, however, beyond this huge central peak in 
ImE(cj), one observes the formation of two almost semi- 
circular structures at uj — ±^ of width ~ D, which can 
be interpreted as an hallmark of the incoherent nature 
of the electronic excitation of the Hubbard band in the 
paramagnetic phase. 



C. Accuracy of the approximation 

After having discussed the results obtained by approxi- 
mating the Hubbard DMFT self-energy starting from the 
t—J one, and having noted that for U ^> D its structure is 
coincident with that of the exact Hubbard model within 
small corrections (roughly of order jjw), here we will ex- 
plicitly compare DMFT data for the optical conductivity 
in order to show the accuracy of our t—J approximation 
for the quantity we are more interested in this paper. 

We report in Fig. [9] our numerical data for the optical 
conductivity of the AF-phase, namely for U = 5D at T = 
calculated inserting the DMFT Hubbard self-energy 
and the correspondent t—J approximated one in Eq. P]l. 
At T = 0, in fact, the possibility to use the Lanczos as 
an impurity solver, allows for a better resolution of the 
multiple peak structure of the Hubbard self-energy and, 
therefore, for a better test case to evaluate the accuracy 
of the t—J approximation for the optical conductivity of 
the Hubbard model. 

In fact, while Lanczos already allows a larger number of 
bath sites than the full exact diagonalization, in the AF 
case for T <C T/v the computational cost of increasing the 
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number of bath sites can be further reduced. Specifically 
it is reasonable to assume that, in the AF ground state, 
also the bath electrons will be almost fully polarized. 

Bearing in mind that the total number of electrons 
with spin up (iVf ) and with spin down (JVj.) are 'good' 
quantum numbers, the Hamiltonian can be decomposed 
in smaller blocks, each associated with a spin sector, 
whose size is given by: 



dira(7iN f ,N l ) = dim(%jv t ) x dim^jvj 
'JV„\ ( N s 



N s -N t 



(12) 
(13) 



where N s represents the total number of electron levels, 
and dim^jv^) the dimension of the Hilbcrt space as- 
sociated with N a electrons of spin a. Since the bath is 
strongly polarized, in the Lanczos algorithm we only need 
to take into account those spin sectors in which the big 
majority of the fermions of the bath has either spin up or 
spin down. That means that we only need to diagonalize 
the smallest blocks composing the Hamiltonian. 

This is very different from the PM phase, where we 



need to take into account mainly those sectors where the 
number of spin up fermions and spin down fermions is 
comparable, and hence the size of the matrices to diago- 
nalize is much bigger. In practice, exploiting the property 
of a fully polarized bath, in the AF phase we can reach 
with a reasonable computational effort a number of elec- 
trons as high as N s = 25, while in the PM case we could 
reach only N s = 14 

Going back to the analysis of Fig. [7J our data show 
that, at U = 5D, the overall agreement between the 
peaky structure of the optical conductivity of the t — J 
approximation and that of the Hubbard model is rather 
good, as both the peak positions and intensity are well 
captured also in the t— J scheme. This result justifies, 
for these values of U, the usage of the DMFT solution of 
the one-hole problem in the t — J model as an accurate 
approximation, e.g., also when considering the finite T 
regime. In fact, when increasing T, Lanczos or full ED 
solution of DMFT are no longer accurate enough for the 
computation of multi-peak spectral functions with suffi- 
cient resolution. 
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For instance, at U — 5D, the maximal magnetization (at 
T = 0) for the Hubbard model is m = 0.979, which corre- 
spond to the t — J solution at /3D — 46.6. This slight differ- 
ence in temperature is naturally reduced with increasing 
U. 

In other words, this means that if the chemical potential 



were rigidly shifted to the value of the first (low-energy) AF 
magnetic peak, this would behave as a coherent quasiparti- 
cle excitation (though strongly renormalized by Z ~ -2j-). 
A similar situation would occur, "mutatis mutandis", in 
the case such an excitation were activated by an exter- 
nal perturbation of frequency oj of the size of the spectral 
gap, e.g., by an electromagnetic field, when measuring the 
optical conductivity a(uj). 



